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ABSTRACT 

We measure the duty cycles for an existing sample of well observed, nearby dwarf 
novae using data from AAVSO, and present a quantitative empirical relation between 
the duty cycle of dwarf novae outbursts and the X-ray luminosity of the system in 
quiescence. We have found that log DC = 0.63(d=0.21) x (log Lx (ergs -1 ) — 31.3) — 
0.95(±0.1), where DC stands for duty cycle. We note that there is intrinsic scat¬ 
ter in this relation greater than what is expected from purely statistical errors . Us¬ 
ing the dwarf n ova X -ray luminosity functions from iPretorius &; Kniggel (|2Q12i) and 
iBvckling et al.1 (|2010fh we compare this relation to the number of dwarf novae in 
the Galactic Bulge Survey which were identified through optical outbursts during 
an 8-day long monitoring campaign. We find a specific frequency of X-ray bright 
{Lx > 10 31 ergs -1 ) Cataclysmic Variables undergoing Dwarf Novae outbursts in the 
direction of the Galactic Bulge of 6.6±4.7 x 10 -5 Mq 1 . Such a specific frequency would 
give a Solar neighborhood space density of long period CVs of p = 5.6±3.9x 10 -6 pc -3 . 
We advocate the use of specific frequency in future work, given that projects like LSST 
will detect DNe well outside the distance range over which p « const. 

Key words: Stars: Novae, Cataclysmic Variables - Stars: Dwarf Novae 


1 INTRODUCTION 

Cataclysmic variables (CVs) are binary star systems in 
which a white dwarf (WD) primary accretes material from a 
low mass secondary star through Roche-lobe overflow. In the 
absence of strong magnetic fields of the WD, the infalling 
material forms an accretion disk, through which angular mo¬ 
mentum can be transfered between particles allowing accre¬ 
tion onto the surface of the WD. 

The material in the disk heats as it builds up, eventu¬ 
ally reaching the ionization temperature of hydrogen. At this 
point, the gas in the disk becomes much more viscous dSmakl 


11994 lL asotal fiooll ). The heat energy from the increase in 
viscosity causes a heating wave to propagate through the 
disk, ionizing the remainder of it. Mass then passes through 
the high viscosity disk at a much higher rate, dumping ma¬ 
terial onto the WD. A sudden brightening of several mag¬ 
nitudes and subsequent decline to the original b rightness 
is observed over th e co urse of days to weeks (Sma ki ll 97ll : 
[Warner k, Natherl ll97ll ). These Dwarf Novae (DNe) occur 
with a frequency tied to the recent mass accretion rate Mi 
and mass transfer rate M 2 of the system, as the disk must 
refill in between outbursts. High mass transfer rate CVs 
(nova-like objects) do not undergo DNe as the disks are con- 
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stantly above the ionization temperature of hydrogen. Low 
mass transfer rate objects such as WZ Sge may go decades 
between recurrences of DN outbursts. 

The mass accretion rate in CVs is also known to be 
tied to X-ray luminosity up to a critical rate at which 
point the Boundary Layer (B L) on the surface of th e WD 
becomes opaque to X-rays d Patters on fc Ravmondl 1 19851 : 
[Wheatley et alJl2000l : lBaskill et al.ll200ll l. Indeed, the X-ray 
emission during a DN outburst can rise at the start of the 
outburst as for SS Cyg until the BL becomes opaque at 
which point the X-rays are quenched. Some DNe never reach 
this critical accretion rate and the X-ray emission remains 
unquenche d through th e duration of the outburst (e.g. U 
Gem; ISwank et a.l.lfl978h . 

The empirical population of CVs can constrain evolu¬ 
tion models; a crucial measurement to make with this end 
in mind is the space density of CVs, p. Models of binary 
evolution predict space densities ranging ov er an order of 
magn itude or more, from 10 -5 — 10 -4 pc -3 (Ide Koollll992l : 
Kolb 1993; Politano 1996). Observational campaigns, largely 
using X-ray selected samples, have made estim ates o f the 


1998|; Schreiber & Gansickc 2003j; 

Auneweroiwit et al. 2006; 

Pretorius et al. 200?!; Rogel et al.l 

20081; 

Pretorius & Knigge 


Pretorius fc Knigge 


d2012ll for long-period CVs was 2.1^ 3 'g x 10 -6 pc -3 which 
assumed that all systems with unknown periods had peri¬ 
ods below the period gap. Studies of the space density of 
CVs have focused on the solar neighborhood, for the simple 
and expedient reason that a complete census of faint CVs 
will be impossible to acquire at large distances. The outlier 
among estimates of space density from observational cam¬ 
paigns, finding a much smaller space density for the dwarf 
nova population (< 5 x 10 -7 pc -3 ), used the OGLE survey 
towards the Galactic Bulge to find dwarf nova outbursts 
llCieslinski et al.ll2003l h This estimate was based on iden¬ 
tifying distant DN outbursts as part of a larg e r wid e-held 
survey as is this work, where ICieslinski et al. (l2003h uses 
the OGLE-II survey and this work is based on the Chandra 
Galactic Bulge Survey (GBS). 

The GBS is an X-ray survey between —3° < l < 3° and 
1° < |b| < 2° with shallow coverage of 2ks (Ijonker et al.l 
l201ll . 120141 ') . The primary goal of the GBS is to conduct a 
census of Low Mass X-ray Binaries in quiescence (qLMXB), 
in order to constrain models of binary evolution and to iden¬ 
tify new systems for which mass determinations are possi¬ 
ble, including eclipsing systems. Towards the goal of classi¬ 
fying each X-ray source in the GBS, multiwavelength pho¬ 
tometric and spectroscopic followup campaigns have been 
undertaken. Some of these focus on classifying new inter- 
esting sources with detailed followup, such as new CVs 


dBritt et al.l l2jll^jRatti et al 


120131') . a possible new black 


hole binary (IBritt et al.ll2013l f. a Slowly Pulsating B star 


( Johnson et al. 2014 

), and 

Binary (Hvnes et al 

20lJ) 


_ _ 2014) ■ In t erms of more comprehen¬ 

sive followup. IBritt et al.l (l2014l 'l focuses on variable opti¬ 
cal counterparts to G BS sources in c ludin g the identification 
of DNe in outburst: Greiss_e£al. (12014!) identifies the in¬ 
frared counterparts: iHvnes et al.l ( 20121 1 cross matches GBS 
sources w ith brigh t star catalogs such as Tycho and ASAS; 
iMaccarone et all (12012! ') cross matches GBS sources with 


existing catalog s of radio sources; iTorres et al.l (|2014l l and 
IWu et all (120151) provide optical spectroscopy of stars with 
an Ha excess near GBS sources. The initial X-ray survey, 
which covered 3/4 of the planned survey area, discovered 
1216 unique X-ray sources with 3 or more detected photons 
in a 2 ks observation. The X-ray observations were kept in¬ 
tentionally shallow in order to keep the number of active 
stars and CVs in the survey low while still reaching the Bulge 
for neutron star and long period black hole qLMXBs. As a 
result, the non-magnetic CVs in the GBS are generally closer 
than Bulge distance. For example, the X-ray bright system 
Z Cam would be detected up to ~ 4500 pc away. 

Our use of an optical variability survey as part of the 
GBS to constrain the space density of CVs in the direc¬ 
tion of the Galactic Bulge mirrors i n some ways that under¬ 
taken with the OGLE-II catalog in ICieslinski et aD d2003lL 
Our survey differs in several key respects. First, our analy¬ 
sis takes into account the likelihood of detecting dwarf no¬ 
vae with recurrence times much longer than the baseline 
of observations. Not allowing for dwarf novae with recur¬ 
rence periods lon ger than the OGL E tim e baseline is the 
likely reason why ICieslinski et al.l (|2003T ) obtained a final 
DN space density much lower than that obtained by other 
methods, or that we find in this paper. Second, the limit¬ 
ing magnitude in our survey is slightly deeper than that of 
OGLE-II, r' = 23 versus I = 22, though this is equalized 
by reddening. Third, our baseline of observations is dramat¬ 
ically shorter than in OGLE-II, with an 8 day baseline com¬ 
pared to OGLE-II’s 4 years of observations, though both 
are shorter than the recurrence time of many CVs undergo¬ 
ing DN outbursts and are therefore far from complete sam¬ 
ples. Fourth, our optically survey ed area is c ur rent ly some¬ 
what smaller, 9 square degrees (IBritt et al.l 120141 ') instead 
of OGLE-II’s 12. Lastly, our optical detections of DNe, at 
the present time, are premised on an X-ray detection of the 
same CV. Recent progress in determining the X-ray luminos¬ 
ity function of CVs combined with a more complete sample 
of well observed, nearby DNe allows our sample to be lever¬ 
aged in ways that were not ava ilable at the time of OGLE-II 
as well. ICieslinski et al.l (12003 1 place an upper limit on the 
space density of CVs towards the Bulge of 5 x 10” ' pc -3 . 
In this work, we reconcile recent estimates of the X-ray se¬ 
lected CV space density with an optical-outburst selected 
dwarf nova space density. 

2 DUTY CYCLE VS X-RAY LUMINOSITY 

2.1 Sample Selection and Measurements for Duty 
Cycle versus Luminosity Relation 

As a basis for determining the relationship between X-ray 
luminosities of CVs and other key p roperties, w e use t he sur¬ 
vey of DN within ~ 200 pc in IBvckling et all ll20ld h This 
sample, although small (13 CVs), includes only objects with 
acc urate distance m easur ements from trigonometric paral¬ 
lax. IBvckling et al . I (i20ldf demonstrate that this sample is 
likely suffering from selection effects against low luminosi¬ 
ties (Lx < 3 x 10 3 °ergs -1 in the [2 — 10] keV band) and 
duty cycles. We therefore limit our fit to the CVs with 
an X-ray luminosity above this cutoff. This sample is not 
distance complete, as it is limited to sources with paral¬ 
lax measurements which were mostly done by the same 
group with Northern hemisphere telescopes dThorstensenl 
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|2003|; iThorstensen et al.l [ 20081; lHarrison et al.| 12004^1 . As an 
exercise, we have fit the entire sample as well and find the 
resulting relation to be consistent with the more narrowly 
defined sample._ 

Since the publication of lBvckling et all (l20ld 'l. the dis¬ 
tance to SS Cyg has been revised downward by t he im¬ 
proved parallax measure in IMiller-Jones et al.l J20131 - ). from 
159 ± 12 pc to 114 ± 2 pc. This results in a drop in the X- 
ray lumi nosity b y a factor of a bout 2 compared to the value 
listed in IBvckling et all d201Cf) . 

Each of these CVs is well observed, with decades of pub¬ 
licly available data in multiple surveys. Where available, we 
use published values of duty cycles extant in the literature. 
For other objects we use AAVSO data in all available filters 
to measure the duty cycle of each CV in this sample, with 
results shown in Table|T] The AAVSO data are generally ac- 
curate to with in 0.2 magnitudes for individual observations 
jHendenll2014l 'l and is much better for many simultaneous 
observations. We define the duty cycle as the number of days 
during which the source is in outburst divided by the to¬ 
tal baseline of observations. We define outburst as an event 
reaching at least 2 magnitudes greater than the quiescent 
magnitude beginning and ending when the object returns to 
0.2 magnitudes above the quiescent level prior to outburst. 
We do not include the standstills in Z Cam’s high state as 
being in outburst; we consider only those times when Z Cam 
is in the low state and undergoing dwarf nova outbursts. We 
only include times at which the sources are visible and being 
observed consistently enough to be sure that no outbursts 
have been missed. Lightcurves of the data ranges used to 
measure duty cycle are presented in Figure |T] 


2.2 Fitting Methods 

Fitting a rel ation b etween the X-ray luminosities in the 
sample from iBvckling et al.l (|201Cll ') and the corresponding 
duty cycles is complicated by the presence of intrinsic scat¬ 
ter around the relation that is not represented by the sta¬ 
tistical errors. Methods for dealing with these difficulties 
al ong with their s t rengt hs an d shor t comin gs are discussed 
in iTremaine et all ll2002l ) and iKelfvl (1200'it l. To account for 
the intrinsic scatter of the data, we add a systematic un¬ 
certainty to all duty cycle measurements that is calculated 
such that fi 2 /v = 1 for the initial fi t by using the MPFI- 
TE XY routine ([ William s et al.|[20ld l in the MPFIT pack¬ 
age (lMarkwardtll200T l . Without using this intrinsic scatter, 
fv = 15 on the final best fit, which is a decidedly poor fit. 
Because the errors in the measurement of X-ray luminos¬ 
ity are somewhat asymmetric, we symmetrize them in log 
space. The changes in the errors this causes are on the order 
of 10%. We use a bootstrapping method to estimate the er¬ 
rors in the fit of the full data set, as the distribution of the 
bootstra pped fit results m atches the error distribution of the 
best fit (iPress et al.l 120071 4. For the bootstrapped fits with¬ 
out including an intrinsic scatter term, l v = H-8 ± 8.4, 
still well above the expected range of X 1 l v = 1 ± ^JFJv. In 
the bootstrapping analysis, the intrinsic scatter in the log of 
the duty cycle is on average amt = 0.31 ± 0.12, which gives 
X 2 /^ = 1.20, which is within x 2 /^ = 1 ± \[CJv for v = 7. 

In order to estimate the errors in the relation, we fit the 
bootstrapped distribution of fits with a Gaussian and use 


Duty Cycle os a Function of X-ray Luminosity 



Figure 2. The best fit for Duty Cycle as a function of X-ray 
luminosity is plotted as a solid line. Dotted lines are the 1 <t errors 
on this relation. 


the width of the Gaussian to estimate the error. The slope 
and intercept have been fit together in each iteration, with 
the intercept set to near the middle of the data to remove 
correlation between the intercept and the slope during the 
fitting. 

The resulting best fit is: 

log DC = 0.63±0.21x (log L x (ergs -1 )-31.3)-0.95(±0.1) 

shown in Figur e [2] Fitting all of the DN in the 
IBvckling et al.l ( 2010 ') sample gives a relation that is con¬ 
sistent with the one above, but which may be subject to 
selection effects against low duty cycle DN with Lx < 
3 X 10 30 ergs -1 , and is given by: 

log DC = 0.65 ±0.08 x (log Lx (ergs -1 ) — 30.5) — 1.50 ±0.14 

This relation is unaffected by the inclusion or removal 
of ASAS J002511±1217.2 (ASAS J0025 henceforth), which 
has only 1 observed outburst, in the sample. The fits are 
shown in Figure [31 

3 THE GALACTIC BULGE SURVEY SAMPLE 

OF DWARF NOVAE 

In the search for variability in candidate optical counterparts 
for the northern 3/4 of th e GBS survey area, described in 
detail in lBritt et alJ ( 2014 1. 8 X-ray sources were observed 
to undergo a DN outburst. Of these 8, only the brightest 
(CX18) has so far been the target of both follow up spec¬ 
troscopic observations and optical monitoring, where it was 
observed to undergo a subsequent DN outburst and to have 
an optical spectru m confirming that it is a non-magnetic CV 
dBritt et al.ll20l.‘j ). 2 others, CX39 and CX87, have been tar¬ 
geted f or spectroscopic fo llow-up which confirms their CV 
nature d Torres el al.llioTil ). The remaining 5 DNe (CX: 81, 
298, 426, 476, 860) have X-ray, infrared, and optical colors 
and flux ratios consistent, where available, with CVs under¬ 
going DN outbursts in addition to the presence of a rise or 
decline of > 1 magnitudes over a few days. Half of the ob¬ 
served DN occur in the brightest 100 X-ray sources out of 
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Table 1. The measured duty cycles and extant X-ray luminosities for nearby dwarf novae. To fit a relation between duty cycle and X-ray 
luminosity, we use only those DN with Lx > 3 x 10 30 ergs -1 . 


System 

log Lx [2-10 keV] 
(ergs -1 ) 

Log Duty Cycle 

Porb 

(hr) 

References* 

BZ Uma 

31.17 ±0.23 

-1.39 ±0.08 

1.63 

i 

HT Cas 

30.79 ±0.21 

-1.85 ± 0.12 

1.77 

i 

SS Aur 

30.982 ± 0.095 

-0.61 ± 0.07 

4.39 

i 

SW Uma 

30.69 ±0.15 

-1.10 ±0.14 

1.36 

i 

U Gem 

30.919 ±0.053 

-0.90 ±0.02 

4.25 

i 

T Leo 

30.81 ±0.13 

-1.17 ±0.16 

1.42 

i 

V893 Sco 

31.70 ±0.31 

-0.74 ±0.12 

1.82 

i 

SS Cyg 

31.888 ± 0.066** 

-0.53 ±0.07 

6.603 

2,3 

Z Cam 

31.79 ±0.33 

-0.38 ± 0.05 

6.98 

4 

GW Lib 

28.7 ± 0.44 

-2.34 ±0.19 

1.28 

1,5 

WZ Sge 

29.85 ±0.12 

-1.85 ± 0.09 

1.36 

1 

VY Aqr 

30.11 ±0.21 

-2.17 ±0.13 

1.51 

1 

ASAS J0025 

30.20 ±0.42 

< -1.92 ±0.15 

1.37 

1,6 


* References: (T)| Hendenl (120141 ) ( 2) lHamilton et al.1 {2007j) (3) I Miller-Jones et alJ (|2013h 
(4) lOppenheimer et al.l~(ll998l) (5 ) lMaza_j£_G onzalcz (1983) IPoimanskil dl997h 
** X-ray luminosities come from iBvckling et al.l (120101) in all cases but SS Cyg 
where we use a revised distance based on radio parallax. 
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Figure 1. AAVSO lightcurves for DNe with parallax measurements within 200 pc. Here we include only those sources with 
Lx > 3 x 10 30 ergs —1 to avoid biases due to data incompletenesses at low Lx- Lightcurves for SS Cyg, U Gem, and Z 
Cam have been binned to 1 day averages for clarity. 
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Duty Cycle os a Function of X-ray Luminosity 



Figure 3. The best fit (solid line) and lcr errors (dotted lines) 
betwe en the duty cycles and X-ray luminosities of the entire sam¬ 
ple in lBvckling et ahl j201C)h ■ The fit to the entire sample agrees 
well with the fit to only the X-ray bright sample, but it is likely 
that there is a strong sampling bias against low duty cycle DN 
below Lx = 3 X 10 30 ergs -1 . 


th e 121 6 source sur vey population. The variability survey 
in lBritt et al.l (120141 ') only examines stars coincident with X- 
ray positions in the GBS catalog; the data exist to study all 
objects in the field, but the full analysis is computationally 
intensive and beyond the scope of this paper. The X-ray de¬ 
tections for the GBS cover the energy range [0.5 — 10] keV, 
and include detections which have at least 3 photons. 

The optical observations were made with the Mosaic-II 
instrument on the Blanco 4m telescope at CTIO in July, 
2012. A description o f the data r educt ion and photome¬ 
try methods is given in lBritt et al.l d2014l ). The observations 
reach an optical depth of r' — 23 and span an 8 day base¬ 
line, which is long enough to notice either the sharp rise or 
slower decay of a DN outburst. 

In addition to t he 8 objects classified as DN outbursts 
in lBritt et al.l |i 20141 ). there is one additional object, CX982, 
which could be a DN outburst with eclipses of the disk by 
the donor star during the outburst. It could also be some 
extreme coronal flaring, so we exclude it from our study 
here. If follow up observations reveal that it is indeed a DN, 
the space density estimated will be too low by ~ 12.5%. To 
account for the chance that this object should be included 
in the sample, and for the smaller but non-zero chance that 
1 or more of our DNe in the sample are, in fact, not DNe, 
we include an additional systematic error of 12% in our final 
estimation of space density. 

4 ESTIMATING SPACE DENSITY WITH 

MONTE CARLO SIMULATIONS 

We can constrain the space density of CVs undergoing DN 
outbursts in the direction of the Galactic Bulge by modeling 
a CV population which f ollows the X-ray luminosity function 
of IPretorius fc Kniggd (120121 ') and the duty cycle relation 
outlined above with various space densities, and comparing 
the resulting population with that seen by the GBS. We 
make several assumptions in this process. 

• CVs trace stellar density in the disk, described by an 


exponential decay in both radius from the center and height 
above the plane of the disk. We have considered 3 possible 
scale lengths of the disk in our simulations, 2.2, 3.0, 3.8 kpc 
to account for the c ontinuing uncertainty in the detai ls 
of Galactic structure djuric et al.l 120081 : iRevle et al.ll2009l ). 
Changes to the scale length did not impact our results in a 
significant way. In addition, we have considered a range of 
scale heights of the disk. Because the scale height assumed 
had some impact on the results, we make the uncertainty of 
the true CV distribution a part of the simulation by adding 
a random component with a variance o f 15 parsecs to a disc 
scale height centered on 120 parsecs ([Pretorius fe Kniggei 
120121) . We assumed a distance to the Galactic Center of 
7620 pc dEisenhauer et al.ll2005l ). 

• The peak absolute magnitude of a DN outburst is usu¬ 
ally between 3 and 6, depending upon the orbital period 
and inclination d War nerlll 9871 : |Pattersonll201ll ). An absolute 
outburst magnitude of 4.5 equ ates to an orbital period of 
4.3 hours using the relation in IPattersonl (|201ll ). We use a 
Gaussian distribution of absolute magnitudes peaking at 4.5 
with a variance of 1.5 magnitudes. To better estimate the 
effects of different distrubutions of absolute magnitudes, we 
reran the simulations with distributions peaking at 5.5 and 
3.5 with the same variance. The changes in the detected 
population of DN outbursts are included in the errors in the 
final space density and specific frequency estimates. 

Because of extinction towards the Bulge, our optical sur¬ 
vey would not recover an outburst with Mv = 4.5 past a 
distance of ~ 5400 pc. Using a range of peak magnitudes as 
described above, we expect « 20% of DNe to be further than 
this, quickly falling off with almost none as far as 8 kpc. At 
5 — 6 kpc in this longitude, the Bulge has started but con¬ 
tributes only ~ 15% of the stellar population in this line of 
sight, though it rises quickly beyond this point. Using mod¬ 
els of stellar density fro m the Bulge popula tion compared 
to the disk population JPasetto et al.l 120121 ). we find that 
the effect of including the Bulge in our models will change 
the number of detected DNe by < 10%, much smaller than 
uncertainties from small number statistics. With a better 
sample of DNe to work with, it would be worth using a 
more accurate model of the Milky Way. We have included a 
10% additional error in the estimate of specific frequency at 
the end of the process to address uncertainty from using a 
simplistic model of the Milky Way. 

• The X-ray spectra of CVs are described by a thermal 
Bremsstrahl ung model with a characteristic temperature of 
15 keV dMukai et al.ll2003l ). Changing the characteristic tem¬ 
perature of the bremsstrahlung model by, e.g., a factor of 2 
up or down makes little difference in the Chandra energy 
band, and thus has little effect on our calculations. 

• Extinction is roughly linear with distance (which is cer¬ 
tainly not really the case since the lines of sight cross spiral 
arms, but this is a small source of error compared with the 
small number sta tistics), and is c alcula ted assuming the ex¬ 
tinction law from lCardelli et al.l d 198 91) . the Vista Var iables 
in Via Lactea reddening maps in I Gonzalez et ahl d2012l) . and 
the relation Nh = 0.58 x 10 22 cm -2 x E(B — V) found in 
iBohlin et al.l d 19781 ). including the observational errors in 
the extinction maps which can be large in these regions. 

• X-ray detections of DNe in the GBS are made outside 
of outburst. This is clearly not accurate - for a “worst case 
scenario”, we would expect for CVs in the GBS like SS Cyg, 
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which is ~ 10 x brighter briefly at the start and end of a DN 
outburst and a high duty cycle (see Table [l]), that about 1/2 
would be detected in outburst due to the larger volume sam¬ 
pled, even though t he brighter luminosity lasts only ~ 10% 
of the duty cycle (iMcGowan et al.l 1200411 . Similarly, more 
systems like U Gem would be detected in outburst than 
in quiescence because the higher X-ray flux lasts through¬ 
out the duration of the outbursts. A survey that did not 
premise the detection of the outburst on an X-ray detection 
would escape this trouble, which could result in somewhat 
overestimating the specific frequency of CVs. 

First, we throw a randomly generated population of CVs 
onto the region of the sky covered by the GBS and calcu¬ 
late the X-ray flux and optical magnitude of a DN outburst 
without correct i ng for extinction. The X-ray luminosities in 
iBvckling et al.l l)2Q10i l are in a somewhat different energy 
band than the GBS data. We use WebPimms to correct for 
this difference using the assumptions about extinction and 
spectral shape listed above. We similarly correct for the dif- 
feren t band used to calcu l ate th e X-ray luminosity function 
from lPretorius fe Kniggd ll2012lf . If the flux values are above 
the GBS detection limits, we calculate the duty cycle based 
on the X-ray luminosity relation above, including the intrin¬ 
sic scatter around this relation seen in the 9 DNe used for 
the fit. At every point that a relation is invoked, the un¬ 
certainty is included as a gaussian random number with a 
variance equal to the uncertainty in the relation. We then 
randomly place an observation at some point in the outburst 
cycle, and if the observation overlaps with a DN outburst, 
the DN is provisionally detected. At this point, we look up 
reddening values for the line of sight each provisionally de¬ 
tected DN is on and calculate extinction in both optical 
and X-ray, including the observational error in the extinc¬ 
tion measures. If the CV remains above the flux limits of 
the GBS at this point, it is detected. 500 such simulations 
are run for different values of space density. We find that 
to generate a population of CVs undergoing DN outbursts 
consistent with those seen in the GBS, we require a space 
density of p = 5.6(±3.9) x 10 -6 pc -3 in the solar neigh¬ 
borhood, scaled up with stellar density as we move closer 
to the Bulge. This is not a measurement of the local space 
density of high duty cycle CVs; it is a measurement of the 
density of high duty cycle CVs several kiloparsecs towards 
the Galactic Bulge, and is expressed in terms of the local 
space density under the assumption that the CV density 
through the Galaxy traces the stellar population. We write 
it in this way to allow a more direct comparison to published 
values of space density. An equivalent, but perhaps easier to 
vizualize, way of writing it is to use the local stellar mass 
density, 0.085 Mq pc -3 jMcMillanl boilll . to state the av¬ 
erage number of CVs per unit stellar mass in the parts of 
the Galactic disk being probed by the optical GBS survey, 
6.6 ±4.7 x 10 -5 Mq 1 . 

This measurement is dominated by the high end of the 
X-ray luminosity function, as lower luminosity CVs undergo 
DN outbursts less frequently and are correspondingly less 
likely to be observed, in addition to the smaller volume 
reached. In particular, in our simulations of the DN popu¬ 
lation in the GBS which best match GBS observations only 
15.2% of detected DN have X-ray luminos ities below the 
Lx = 3 x 10 3O ergs _1 completeness limit of IBvckling et all 


d2010h . as shown in Figure []] which would correspond to 
roughly 1 of the 8 detected DN outbursts in the GBS. As 
such, the exact shape of the assumed luminosity function at 
low to moderate luminosities has little effect on the resulting 
population, and the primary uncertainty comes ultimately 
from the small number of systems involved at both ends of 
the simu lation. For e xamp le, using the luminosity function 
from IBvckling et al.l (|2010l ) instead does not have a signifi¬ 
cant impact on the resulting specific frequency. We also have 
very little in the way of a lever arm to constrain the relative 
size of the X-ray faint DN population to the X-ray bright 
DN population. The space density we measure is a factor 
of 2 higher than, but is still co nsistent with, the measure - 
ment of local long period DN in lPretorius fe Kniggel ll2012h . 
2.lj; 3 ; 3 x 10 -6 pc -3 . 

5 DISCUSSION 

Studies of the space density of CVs have concentrated on 
the local population, so that the changing density of the 
stellar environment as one moves outward radially in the 
disk is unimportant and generally not considered. In the 
GBS, however, we likely observe CVs undergoing DN out¬ 
bursts out to 4 — 5 kpc. The radial and vertical exponential 
decay of the stellar density is therefore very important in 
our estimate of space density in this line of sight. GBS fields 
are also particularly close to the plane compared to other 
measures of space density of CVs which have tried to avoid 
extinction as much as possible by looking out of the plane 
of the Milky Way. The concept of a single “space density” 
for CVs is one that only makes sense in an extremely local 
environment where the density is relatively uniform, and is 
only a useful concept for the population at large in terms of 
a scale factor to the stellar density. 

While the non-uniform structure of the Galaxy is taken 
into account when calculating the space density of CVs, a 
more straightforward and ultimately more useful concept 
may be the specific frequency of CVs, which offers a more 
direct comparison to population models throughout differ¬ 
ent regions of the galaxy than the solar environment. The 
optical photometric follow-up to the GBS jBritt et al.ll2014l l 
reached a depth of r' = 23. The optical depth of the survey 
combined with the bright X-ray luminosities expected from 
the very fact of the outburst detection in such a short base¬ 
line results in the distribution shown in the right panel of 
Figure]]]] This work is therefore a probe of the density of non¬ 
magnetic CVs at distances from Earth much greater than 
any previous survey. Comparing our res ult to the local space 
densit y of long period CVs reached in IPretorius fc Kniggel 
(l2012h . there is some hint at a higher specific frequency of 
CVs as one moves in towards the center of the Galaxy; this 
is not yet a statistically significant effect, but we have a 
small number of DNe with which to work and make no at¬ 
tempt to factor in the distances of DNe in the GBS as a 
result. There is both a theoretical and observational basis 
to expect a ch ange in specific frequ ency in different Galactic 
environments. iKundu et al.l l|20(M ) show that an increase in 
the metallicity of Globular Clusters results in an increase in 
the specific frequency of Low Mass X-ray Binaries in those 
clusters. A possible reason for the increase is that donor 
stars with low metallicities do not develop convective zones 
at the same stellar mass, and so do not make contact un¬ 
til shorter orbital periods have been reached through other 
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Figure 4. Left: A histogram of the X-ray luminosities of the simulated population of DNe detected by the G BS u sing the relation between 
duty cycle and X-ray luminosity derived in Section 12.21 Systems below the break in the iBvckling et al.l |2010| ') luminosity function are 
plotted in red with hatching at —45°, while those above are in blue with hatching at +45°. 84.8% of systems are above the break. The 
GBS sample of DNe, therefore, likely is composed primarily of systems of comparable luminosity to those used to fit the relation between 
duty cycle and X-ray luminosity. Also and importantly, the sample of DNe detected in the GBS is relatively insensitive to the precise 
shape of the luminosity function, with space density of X-ray luminous CVs playing the major role. This is a result of the short baseline 
of GBS observations (only 8 days); longer surveys such as LSST will detect non-magnetic CVs with much lower duty cycles, offering 
a window into the distribution of the X-ray faint CV population. Right: The simulated distribution of distances of detected DNe. The 
optical depth of the survey combined with the bright X-ray luminosities expected from the very fact of the outburst detection in such 
a short baseline results in a distribution that falls off after ~ 6kpc. The primary limit to how far we can detect DN outbursts is the 
optical brightness of the outburst compared to our survey depth. 


mea ns of angular mome ntum loss such as magnetic brak¬ 
ing lllvanova et al.ll200a ). We consider it likely, therefore, 
that the specific frequency of CVs should increase towards 
the Galactic Bulge as the metallicity increases, though these 
data are insufficient to show it. 

This present work is severely limited by the small num¬ 
ber of DNe available for the analysis in both the high lumi¬ 
nosity Byckling sample and in the sample of DNe outbursts 
detected in the GBS in both optical and X-ray wavelengths. 
Though the DNe in the GBS sample are all detected in the 
X-ray, this method is primarily based on DN outburst detec¬ 
tion, using the X-ray luminosity function and the connection 
between X-ray luminosity and Duty Cycle to connect to a 
predicted outburst rate. Our agreement with other space 
density estimates for DNe is a step forward for DN outburst 
de tection methods, a s other DN outburst surveys such as 
ICieslinski et al.l |2003l ) are in disagreement with other meth¬ 
ods. Our result does not resolve any apparent tension be¬ 
tween empirical and theoretical estimates of CV space den¬ 
sity discussed in |§ 1| as it agrees well with the majority of 
other empirical methods. 

LSST will offer a much better sampling of DNe sys¬ 
tems. LSST will routinely go deeper in all parts of the sky, 
with frequent sampling and a long baseline that guaran¬ 
tees the detection of thousands of DN outbursts, including 
those with low duty cycles. Indeed, some large survey work 
revealing hu ndreds of new DN e outbursts is already forth¬ 
coming (e.g. iDrake et al.l 120141 ). Through this work, using 
large survey data to measure the frequency and brightness 
of outbursts for thousands of systems can yield an estimate 
of the relative proportion of CVs over all parts of the sky. 
Given the consistency of the optical brightness of DN out¬ 
bursts (iPattersonll201ll ). maps in 3 dimensions to distances 


of a few kpc could be made with no observations beyond 
what LSST is going to take anyway. A similar analysis with 
Pan-STARRS is complicated by the observing cadence on a 
similar or longer timescale than the duration of many DNe 
outbursts, many of which may only last a matter of days. 

The primary uncertainty of our estimate here comes 
from the small number of systems used, both to establish 
a quantitative relationship between duty cycle and X-ray 
luminosity and to match a synthetic population t o our 8 
DNe in the GBS optical campaign of iBritt et al.l (l2014|j . 
The former is difficult to improve at present, as the number 
of X-ray luminous CVs for which highly accurate distances 
can be measured is small. However, the GAIA mission will 
undoubtedly improve this situation dramatically in several 
years’ t i me, roughly concurrently with results from LSST 
(iBarstow et al.ll20l4 b 

Neither the DN duty cycle nor the X-ray luminosity 
in quiescence are entirely reliable tracers of the secular mass 
accretion rate of the system (for which A/i = —M 2 in conser¬ 
vative mass transfer), as CVs at the same orbital period can 
exist in states of greatly differing mass tr ansfer rates (e.g. th e 
nova-like below the per iod gap, BK Lvm lZellem et aLl 2009). 
iBiining &; Ritterl ||2004| ) suggest that this difference could be 
due to feedback cycles from irradiation driven mass transfer 
on the scale of 10 4 ~ 5 years, but regardless of the mecha¬ 
nism, the disparate instantaneous mass transfer and accre- 
tion rates at equal o rbital period is an observational fact. 
IPatterson et all (|2013l ) found that BK Lyn is transitioning 
to lower mass transfer rates, moving from a nova-like CV 
to a Z Cam-like state that alternates between a high and a 
low state, perhaps "calming down" after an historical nova 
in 101 A.D. quickly altered the equilibrium of the system. 
Z Cam itself could be in a similar transitional period af- 
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ter a nova eruption rtShara et alj|2007h . If these periods of 
high mass transfer only last a few thousand years after a 
nova, then the majority of systems should have a strong 
relat ionsh ip betw een mass transfer rate and orbital period. 
IPattersonl (1201 lh found that the recurrence time of DNe be¬ 
low the period gap was strongly a function of mass ratio, 
with T rec = 318(±30) days (g/0.15© 2 ’ 63± ° 17 . The mass ra¬ 
tio also determines the Roche geometry — and is strongly 
related to the WD temperature, itself a function of WD 
mass and accretion rate (see Figure 10 of IPattersonl l201lh . 
There is no theoretical prediction of the quiescent instan¬ 
taneous accretion rate (Mi) based on system parameters, 
though these observational signatures are suggestive. Our 
result may therefore be unsurprising to the community, but 
is the first quantitative relationship between the duty cy¬ 
cle and X-ray Luminosity, and so directly links the average 
mass transfer rate to the instantaneous accretion rate in 
quiescence. There is a strong selection effect on orbital pe¬ 
riod in our sample, as wel l, since long perio d systems have 
more massive donor stars lIFrank et al.|[2002f ) and there is a 
maximum WD mass. Above the period gap, the mass ratio 
q must be higher than q ss 0.25, and will have relatively 
short recurrence times. Long period systems will therefore 
be significantly over-represented in a sample of X-ray lumi¬ 
nous, high duty cycle DNe such as those in the GBS, though 
short-period systems are also sure to be present. 
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